Enhanced mesoscopic fluctuations in the crossover between random matrix ensembles 
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In random-matrix ensembles that interpolate between the three basic ensembles (orthogonal, uni- 
tary, and symplectic), there exist correlations between elements of the same eigenvector and between 
different eigenvectors. We study such correlations, using a remarkable correspondence between the 
interpolating ensembles late in the crossover and a basic ensemble of finite size. In small metal grains 
or semiconductor quantum dots, the correlations between different eigenvectors lead to enhanced 
fluctuations of the electron-electron interaction matrix elements which become parametrically larger 
than the non-universal fluctuations. 
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Random matrix theory has focused on the study of 
three ensembles of Hamiltonians: the Gaussian Uni- 
tary Ensemble (GUE), the Gaussian Orthogonal En- 
semble (GOE), and the Gaussian Symplectic Ensemble 
(GSE). These describe the statistics of single-particle en- 
ergy levels and wavefunctions of disordered metal grains 
or chaotic quantum dots with the corresponding sym- 
metries; GUE if time-reversal symmetry is broken, and 
GOE or GSE if time-reversal symmetry is present and 
spin-rotation symmetry is present or absent, respectively. 
In these three basic ensembles, eigenvector elements are 
Gaussian complex/real/quaternion random numbers; el- 
ements of the same eigenvector and of different eigenvec- 
tors are all statistically independent Q|. 

Disordered or chaotic systems with partially broken 
symmetries show a variety of phenomena that go beyond 
a mere "interpolation" of descriptions based on the GOE, 
GUE, and GSE alone. For example, in a quantum dot, a 
weak magnetic field causes long-range wavefunction cor- 
relations and a non-Gaussian distribution of "level 
velocities" , derivatives of energy levels with respect to, 
e.g., a shape change of the dot ||. Both effects are ab- 
sent without a magnetic field (in the GOE), or when 
the magnetic field is strong enough to fully break time- 
reversal symmetry (in the GUE). In a metal grain, weak 
spin-orbit interaction induces mesoscopic fluctuations of 
the g-tensor [^|J^], which does not fluctuate in either the 
GOE or the GSE. Further, as we'll show below, in a weak 
magnetic field or for weak spin-orbit scattering, matrix 
elements of the electron-electron interaction exhibit fluc- 
tuations that are parametrically larger than in each of 
the three basic ensembles. 

The underlying reason for these phenomena is that 
eigenvector elements are not independent in (random- 
matrix) ensembles that interpolate between the three 
basic symmetry classes: There exist both correlations 
within the same eigenvector and, as we show in 

this letter, between different eigenvectors. To study the 
eigenvector correlations in such crossover ensembles, we 
will make use of a surprising relation between the eigen- 
vector statistics late in the crossover from class A to class 
B and that of finite-sized matrices in class B (where B is 



the class of lower symmetry) . Examples of such a relation 
were known for the statistics of a single eigenvector. For 
example, in the GOE-GUE crossover, which is described 
by the N x N random hermitian matrix (with N taken 
to oo at the end of the calculation) 



H 0V (N,a) = H (N) 



(1) 



the distribution of the "phase rigidity" | u» T w | 2 Q] of a sin- 
gle eigenvector v is the same as in the finite-sized M x M 
GUE ensemble with M — 2c? if a is large. In Eq. (g), 
Ho(N) and Htj(N) are N x N matrices taken from the 
GOE and GUE, respectively, with equal variances for the 
matrix elements. A similar correspondence occurs for the 

§ tensor of a Kramers doublet in the GOE-GSE crossover 
,Q . Our main finding is that such a correspondence ex- 
tends to the correlations between different eigenvectors. 

In this paper we will accomplish four tasks, (i) We 
show numerically that the relation 



ffou(JV,a)Hff D (M), M = 2a 2 



(2) 



between the GOE-GUE crossover Hamiltonian Hon (a) 
for large a and N and a finite-sized MxM GUE Hamilto- 
nian extends to correlations between eigenvectors. Just 
as in critical phenomena, where simple power laws un- 
fold into universal scaling functions as you flow away 
from the critical point, here a rich theory of correlations 
unfolds in the crossover region. We wish to point out 
that this principle applies not only to the GOE-GUE 
crossover, but also, e.g., to the GOE-GSE crossover, or 
to wavefunctions in two coupled quantum dots, which 
are described by a random Hamiltonian interpolating be- 
tween two independent GUE's and one GUE of double 
size ||. (ii) We show that, for large a, the universal- 
ity classes are actually curves in the (l/a,l/N) plane, 
reminiscent of renormalization-group flow trajectories 
fiof . (iii) We calculate correlations between eigenvectors, 
based on the surmise (|J) and diagrammatic perturbation 
theory, (iv) We calculate how the inter-eigenvector cor- 
relations in the crossover region affect matrix elements 
of the electron-electron interaction in a quantum dot or 
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metal grain in a weak magnetic field, and predict a signif- 
icant enhancement of fluctuations compared to the basic 
ensembles. 

Let us now consider the joint distribution P({w M }) of 
n eigenvectors v^, p = 1, . . . , n, for the example of the 
GOE-GUE crossover Hamiltonian (|l|). Throughout the 
entire GOE-GUE crossover, the distribution of the eigen- 
vectors is invariant under orthogonal transformations. As 
a consequence, the joint distribution P({v^}) is com- 
pletely determined by the distribution of the orthogonal 
invariants (2|^] 

Pnv = Pun = v^v u) p,u = l,...,n, (3) 

where the superscript T denotes transposition. Hence 

~ n 

P(M) = J Hd P ^P({p M u}) 

n 

x JJ 6(vlv v - 5^)8{vlv v - Pixv ). (4) 

For the physically relevant case of large N, Eq. (Q) implies 
that the eigenvector elements f Mm , m = 1, . . . , N, have 
a Gaussian distribution with zero mean and 

The subscript (. . .}„ indicates that the average is taken at 
fixed Pliv . For the full ensemble average one has to per- 
form a subsequent average over the p M „ with the distribu- 
tion P{{pnu})- We can find P{{pp,v}) from the surmise 
that, for a > 1 and for eigenvectors whose energies 
are all inside a window of size < a 2 A, A being the level 
spacing of the Hamiltonian H{a), the joint distribution 
of the p^ v is the same as for a GUE Hamiltonian of fi- 
nite size M — 2a 2 . Thus the p^ v are independently and 
Gaussian distributed with zero mean and with variance 

(\p^\ 2 ) = (l + S^)/M, M = 2a 2 . (6) 

Together, Eqs. (0)-(@) fix the joint distribution of eigen- 
vectors in the crossover ensemble close to the GUE. For 
the single-eigenvector distribution, Eqs. (@)-(@) repro- 
duce the a 3> 1 limit of the exact solution of Sommers 
and Iida || . The fact that the phase rigidity | ^ | 2 of a 
single eigenvector is a fluctuating quantity is the prime 
cause of the correlations between elements of one eigen- 
vector ^|,|| ; It is the existence of nonzero and fluctuating 
p^u for p 7^ v that causes the correlations between dif- 
ferent eigenvectors. 

We now proceed to present arguments in support of our 
surmise. We consider the eigenvectors v „ (/i = 1, . . . , n) 
with energies within a distance <C a 2 A from a reference 
energy e rc f, sorting them by increasing energy. We then 
consider how each of these eigenvectors is built up from 
the eigenvectors o v of the unperturbed Hamiltonian Hq . 
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FIG. 1. Left panel: Eigenvalues for one realization of 
Ho\j(a). The shaded region marks the energy window of size 
~ M(q)A = 2a 2 A for which the eigenvalues are kept in the 
effective M x M GUE Hamiltonian. Right panel: (\p^\ 2 ) 
as a function of the distance p — v ~ (e M — e v )/A between 
eigenvalues, for a = 4.0 (dashed line, left panel). Solid curve: 
Eq. (|^). Data points: numerical calculation for N = 400. 

Contributions from eigenvectors o„ with energy e v far 
away from e rc f can be described using perturbation the- 
ory in the crossover parameter a. If the energy difference 
£ref — is sufficiently large, the admixture of o u to any 
of the vectors v^a) of interest is small and can be ne- 
glected. On the other hand, eigenvectors o u with energy 
e v close to £ rc f contribute non-perturbatively for large a. 
Upon increasing a from zero, the eigenvectors tv( a ) m 
this energy range have undergone several avoided cross- 
ings, and the unperturbed eigenvectors o M have roughly 
equal weights in each of the vectors v^(a) in our set. 

It is on this heuristic picture that our surmise for an ef- 
fective description of the eigenvector statistics for large a 
is based: We only retain those eigenvectors of the unper- 
turbed Hamiltonian Ho that are relatively close in energy 
and hence all contribute roughly equally, see Fig. |l| for 
a cartoon. Since the time-reversal symmetry breaking 
perturbation in Eq. ([!]) is strong for these eigenvectors, 
the matrix elements between them form a random hermi- 
tian matrix of the GUE. Denoting the effective number of 
contributing unperturbed eigenvectors as M(a), we thus 
reduce the problem of finding the distribution of the or- 
thogonal invariants p^ v for the N x N crossover Hamil- 
tonian (Q) to that of finding the distribution of the p^ v 
for the much smaller GUE Hamiltonian of size M{a). To 
calculate M(a) in terms of ./V and a, we turn to the exact 
solution for the single-eigenvector distribution obtained 
in Refs. §-§, and find fn| 

M[a) = a 2 N(a 2 + 2N)/{a 2 + N) 2 . (7) 

For large N this simplifies to M(a) = 2a 2 , in agreement 
with Eq. @. 

By our surmise, the distribution of the orthogonal in- 
variants should depend on the effective matrix size M(a) 
only, not on a and N individually, as long as N and a 
are large. We have verified this by numerical calculation 
of the averages {\pp V \ 2 ) for different points along a curve 
of constant M(a) in the (l/N,l/a) plane. The results 
of such a calculation are shown in Fig. for p = v, for 
neighboring eigenvectors (p = v+1) and for next-nearest 
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FIG. 2. Curves of constant effective GUE size M(a), Eq. 
([?]), in the (I/a, 1/N) plane for the N x N crossover Hamil- 
tonian (§). Top to bottom: M = 20, M = 30, M = 50, 
M = 100, and M = 200. The horizontal and vertical axes 
correspond to the pure GUE and to the N — > oo crossover 
Hamiltonian, respectively. Inset: (Ip^l 2 ) for the points indi- 
cated at the M = 50 curve in the main panel. Circles: p = v\ 
Squares: p = v + 1 (eigenvectors with neighboring energy 
levels); diamonds: p = v + 2 (next- nearest neighbors). The 
dashed lines indicate the surmise of Eq. ([]) . 

neighbors (p = v + 2). We have also verified that the dis- 
tribution of the Out, is indeed Gaussian (not shown). 

The surmise (g) is expected to be valid as long as 
only eigenvectors taken from an energy window of width 
<C M{a) A = 2a 2 A are involved. If the energy differences 
between eigenvectors become of order a 2 A or larger, the 
eigenvectors v^(a) do not share the same unperturbed 
eigenvectors o„, and we thus expect that they become un- 
corrected. A quantitative description of eigenvector cor- 
relations at energy separations 3> A can be obtained us- 
ing diagrammatic perturbation theory. The only nonzero 
second moment is (\pp,v\ 2 ), which can be computed from 



\Piu>\ ) = 



A2_ 



E 

Si ,S 2 =± 

T 



SlS 2 



x (trG T (e A1 + is 1 6)G{e v +is 2 S)) 



(8) 



where G{z) — l/(z — Hon), £ is a positive infinitesi- 
mal, and the eigenvectors and v v have energies Ep 
and e„, respectively. Calculating the averages using the 
technique of Ref. p2[, we find, if p ^ v, 



(l/Vl 



2a 2 



4a 4 +7T 2 (£ M -£ I ,) 2 /A 2 



(9) 



A similar result for parametric correlations inside a basic 
random-matrix ensemble was derived in Ref. [ p^[ . The 
right panel of Fig. |l| shows (\pp V | 2 ) as a function of Sp — e v 
and a numerical calculation of the same quantity. 

The GOE-GUE crossover that we considered here de- 
scribes the wavefunction statistics in, e.g., a chaotic 
quantum dot or a disordered metal grain in a weak 
magnetic field. Wavefunction distributions have immedi- 
ate experimental relevance for the spacings, widths, and 



heights of Coulomb blockade peaks in the conductance of 
metal grains or quantum dots flij . Correlations between 
wavefunctions of neighboring energy levels cause corre- 
lations between the heights and widths of neighboring 
conductance peaks. Wavefunction distributions also in- 
fluence the positions of Coulomb blockade peaks through 
their role in the distribution of electron-electron interac- 
tion matrix elements |l5| ] , which we now discuss in detail. 
The interaction matrix element U^pu is defined as 



U, 



= / dr\dr-2\J(f\ — f 2 ) 



Mri)Mr2)Mr2)*MnT, (10) 

where U(r) is the electron-electron interaction poten- 
tial and <fifi(r) the wavefunction for an electron in level 
e„. For example, the difference of interaction matrix el- 



ements U u 



Uptptoo gives the spacing between peak 



positions corresponding to different nonequilibrium con- 
figurations (levels v and o unoccupied, respectively) in 
tunneling spectroscopy of small metal grains ]l6| . 

In a metal grain or quantum dot, the interaction can 
be approximated by an r-independent part and a local 
interaction U loc (f) = XAV5(r), where A is the mean 
level spacing, V is the sample volume, and A is a param- 
eter of order unity governing the strength of the local 
interaction. The spatially constant interaction leads to 
a charging energy and does not show mesoscopic fluc- 
tuations. For the GOE (no magnetic field), the matrix 
elements of U loc have an average given by |14| 

(E#SJ = AA (<W<W + &»P&»o + Spjpa)- (11) 

If time-reversal symmetry is broken by a magnetic field 
(i.e., in the GUE), the last term in Eq. Q is left out 0. 
In both the GOE and GUE, fluctuations of the interac- 
tion matrix elements U x °„ pa and corrections to Eq. ( |Tl| ) 
are nonuniversal and small as (at most) g~ x ^ 2 , where g is 
the dimcnsionless conductance of the sample. Equation 
(|ll| ) can be reproduced from random-matrix theory if the 
wavefunctions (j)^ (r) are replaced by eigenvectors and 
the integration over space is replaced by a summation 
over the vector indices. 

How are the interaction matrix elements distributed 
in the presence of a weak magnetic field? If we are not 
interested in the non-universal (l/<?) corrections, that 
question can be answered using the eigenvector distribu- 
tions for the GOE-GUE crossover that we derived above. 
First, upon increasing the magnetic field, there is a grad- 
ual suppression of the last term in Eq. (|i"T|). However, 
as a result of gauge invariance (wave-functions can be 
multiplied with an arbitrary phase factor), the average 
of the last term in Eq. (|ll]) is formally zero throughout 
the crossover if ji ^ p, and this effect only shows up 
in the fluctuations of U^ pa -. Second, the appearance of 
inter-eigenvector correlations enhances the average of all 
interaction matrix elements. Using Eq. (|^), we find 
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FIG. 3. Root-mean-square fluctuations of the interaction 
matrix element U pvpa for four consecutive levels fi = v — 1 
— p — 2 = a — 3. The dashed line shows the large-Q asymp- 
tote of Eq. (|l4|). The solid line is obtained from numerical 
generation of 400 x 400 GOE-GUE crossover matrices, using 
(\U pvptT \ 2 ) — (AA) 2 (|p M „| 2 |p pCT | 2 ). (Direct numerical calcula- 
tion of Up. vpa suffers from large finite- iV corrections.) 

(U l ™ pa ) = \A(5 pp 5 ua + 5^5 vp + ( P ^ Ppa )). (12) 

For large a, the average (p pv p pa ) can be done with the 
help of Eqs. (^J) and (0), with the result 

X { 1+ 4a± + n 2 (e p -e„) 2 /A 2 )- ( 13 ) 

Third, the inter-eigenvector correlations enhance the fluc- 
tuations of the interaction matrix elements. This is best 
illustrated by the expectation value (\U ^ pa \ 2 ) with all 
four indices p,, p, and a different, 

(|C/XJ 2 > = (AA) 2 (|p^| 2 | PpCT | 2 ) = (AA) 2 /(2a 2 ) 2 . (14) 

The first equality in Eq. ( |l4| ) holds for all a, the sec- 
ond one only if a 3> 1 and the four eigenvalues e M , e v , 
e p , e a are within a distance <C cv 2 A of each other. We 
have numerically calculated (|l/^° c pct | 2 ) for four neighbor- 
ing wavefunctions as a function of a throughout the en- 
tire GOE-GUE crossover, see Fig. ||. 

Note that the enhanced fluctuations appear in ma- 
trix elements that are not commonly associated with the 
presence of time-reversal symmetry. Although the fluc- 
tuations are small if a ^> 1, they can be significantly 
larger than the non-universal fluctuations that vanish as 
g~ 2 [for Eq. (Q)]. These extra large fluctuations can be 
an additional source of fluctuations of Coulomb-blockade 
peak positions, and may dominate over the non- universal 
sources of fluctuations. 

The origin of the eigenvector correlations and the en- 
hanced fluctuations of interaction matrix elements can 
be sought in the existence of the large parameter a 2 that 
plays a role similar to the dimensionless conductance g 
in the pure ensembles. The parameter a 2 can be iden- 
tified as the ratio of the Heisenberg time th = 2ttH/A 
and the time tqu needed to acquire a flux quantum [jlJJ . 
Late in the crossover, GUE physics ranges from the mean 
level spacing A up to the scale fi/rou- In the pure GUE, 



however, validity of random-matrix theory ceases only at 
the higher energy scale S/r erg , where r crg is the ergodic 
time. The role of the large parameter g = tjt / i~ cr g , which 
governs wavefunction correlations and interaction matrix 
element fluctuations in the "pure" GUE and GOE is thus 
played by a 2 ~ th/toh in the GOE-GUE crossover. 
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